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ABSTRACT 


Reflection  and  transmission  of  electromagnetic  plane 
waves  at  the  surface  of  parallel-plate  metal  structures  are 
studied  both  theoretically  and  experimentally.  This  report, 
Part  I  of  the  study,  is  tte  theoretical  presentation;  Part  II 
will  deal  with  apparatus  and  experimental  findings.  The 
study  is  confined  to  rays  with  polarization  parallel  to  the 
edges  of  the  plates  and  a  plane  of  incidence  perpendicular  to 
these  edges. 

The  theory  of  Carlson  and  Heins  is  extended  and  made 
amenable  to  numerical  predictions.  Tables  and  graphs  are 
provided  for  the  coefficients  pertaining  to  an  air -metal -plate 
medium  interface.  A  method  is  developed  for  the  analysis  of 
reflecting  and  transmitting  properties  of  slabs,  or  media  of 
finite  depth. 


PROBLEM  STATUS 

This  is  an  interim  report;  work  on  this  problem  is 
continuing. 
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PHYSICAL  OPTICS  OF  METAL-PLATE  MEDIA; 
PART  1,  THEORETICAL  CONSIDERATIONS 


INTRODUCTION 

This  report  deals  with  certain  aspects  of  wave  propagation  in  the  presence  of  the 
parallel-plate  metal  structures  employed  in  microwave  lens  design.  Tliese  structures 
consist  of  a  set  of  parallel,  equidistant  conducting  plates.  They  are  so  oriented  that  the 
polarization  of  the  electric  vector  of  the  Incident  waves  is  parallel  to  the  plates  and  the 
plate  separation  is  chosen  between  one-half  and  one  wavelength.  When  this  is  the  case 
the  phase  velocity  within  the  metal-plate  structure  is  greater  than  in  free  space.  This 
fact  has  been  utilized  by  Stiitzer ,*>*»''  Kock,*  and  others  for  the  construction  of  microwave 
lenses.  The  lenses  so  constructed  are  unusual  in  their  appearance  because  the  metal- 
plate  medium  has  an  effective  index  of  refraction  less  than  one,  consequently  convex  lenses 
are'dlverging  and  concave  ones  converging. 

The  geometrical  optics  of  metal-plate  media  has  already  been  thoroughly  explored. 

The  present  study  is  directed  toward  the  exploration  of  the  physical-optical  aspects  of 
these  media. 

Articles  and  reports  dealing  with  microwave  lens  design  indicate  that  the  temptation 
to  treat  metal-plate  media  formally  as  ordinary  dielectric  materials  has  been  great.  While 
®ost  relations  of  geometriC&l  optics  can  be  carried  over  formally  from  dielectrics  to 
metal-plate  media,  this  is  not  true  of  the  relations  of  physical  optics.  The  different  be¬ 
havior  of  the  latter  media  is  in  a  large  part  caused  by  the  fact  that  energy  is  stored  at  the 
boundary  surface  of  the  metal-plate  medium  while  no  such  storage  takes  place  at  the  sur- 
lace  of  the  dielectric. 

The  objective  of  this  study  is  to  examine  the  reflecting  and  transmitting  properties  of 
metal-plate  media.  The  ground-work  in  this  field  has  already  been  laid  by  otoer  investiga¬ 
tors.  In  1939  Stiitzer®  obtained  approximate  expressions  for  the  reflection  and  transmission 
coefficients  for  normal  incidence.  The  wave  constructions  which  led  to  these  coefficients, 
however,  do  not  satisfy  the  proper  boundary  conditions.  Several  years  later  Carlson  and 
Heins*  attacked  the  problem  anew.  They  solved  the  field  problem  rigorously  for  any  angle 
of  incidence,  requiring  only  that  the  plane  of  incidence  be  perpendicular  to  the  metal  plates. 
By  the  use  of  the  Fourier  transformation,  Carlson  and  Heins  obtained  analytic  expressions 
or  the  reflection  and  transmission  coefficients  which  are  valid  when  the  incident  radiation 
to  subject  to  a  number  of  restrictions. 

experimental  program  has  been  initiated  at  the  Naval  Research  Laboratory  and  at 
State  College  with  the  purpose  of  testing  the  predictions  of  the  theory  of  Carlson 
^  Heins.  The  experimental  program  and  its  results  are  the  subject  of  a  separate  report. 
®to:sory  examination  of  the  work  of  Carlson  and  Heins  is  sufficient  to  convince  one  that  a 
°*tolderable  amount  of  work  is  required  to  bridge  the  gap  between  the  mathematical  formalism 
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of  their  theory  and  the  variables  accessible  to  the  experimenter.  The  theory  predicts  th* 
reflection  and  the  transmission  of  a  plane  electromagnetic  wave  on  a  single  boundary  ^ 
face  separating  free  space  and  a  metal-plate  medium  extending  (in  depth)  to  infinity/ 
case  of  a  medium  of  finite  depth,  the  reflection  from  the  ‘back”  surface  combines  Mth  ^ 
at  the  front  surface.  An  observer  can  observe  only  the  combination  of  these  reflectiote^ 
unless  he  can  eliminate  all  reflection  at  the  back  surface  by  the  introduction  of  prop^ 
absorbing  materials.  In  an  experimental  test  of  the  theory  it  is  then  necessary  either  to 
compute  reflection  and  transmission  coefficients  for  media  of  finite  depth  from  the  slngi# 
surface  coefficients,  or  it  is  necessary  to  eliminate  back  surface  reflections  by  the  u^^^ 
absorbing  materials.  Both  methods  have  been  employed,  the  first  one  primarily  at  the 
Naval  Research  Laboratory,  the  second  at  Oregon  State  College. 


The  first  method  requires  the  knowledge  of  the  phase  as  well  as  the  magnitude  of  the 
coefficients  associated  with  a  single  surface  transition.  While  phase  relations  between 
incident,  transmitted  and  reflected  waves  are  trivially  simple  for  ordinary  lossless  dlelec. 
trie  materials,  quite  the  opposite  is  true  of  metal-plate  media.  The  expressions  for  phase 
relations  derived  by  Carlson  and  Heins  are  of  the  form  of  infinite  series;  their  numerical 
evaluation  presented  a  number  of  difficulties.  In  view  of  this  fact,  the  technique  of  the  com. 
putation  is  described  in  this  report  and  tables  of  the  coefficients  computed  are  appended. 


Expressions  for  the  reflection  and  transmission  coefficients  were  derived  by  Carlson 
and  Heins  for  the  case  when  to  every  incident  beam  corresponds  a  single  reflected  beam. 
This  condition  limits  the  angle  of  incidence  below  a  critical  value  determined  by  the  ratio 
of  wavelength  to  plate  separation.  When  the  critical  value  is  exceeded,  the  edges  of  the 
plates  act  as  a  grating  and  diffracted  beams  arise.  By  extending  the  theory  of  Carlson 
and  Heins,  expressions  have  been  found  for  the  reflection  and  transmission  coefficients  in 
the  presence  of  a  diffracted  beam.  Tables  and  graphs  have  been  prepared  to  show  the  en¬ 
ergy  distribution  between  various  beams.  The  tables  are  not  as  extensive  as  might  be  de¬ 
sired  because  of  limitation  in  time  and  manpower.  Nevertheless  those  angles  of  incidence 
which  are  of  practical  value  are  covered. 


Having  obtained  the  basic  coefficients  pertaining  to  a  single  air -metal-plate  medium 
interface  it  became  necessary  to  consider  the  problem  of  composing  reflections  originating 
at  the  front  and  back  surfaces  of  a  slab,  or  medium  of  uniform  depth.  This  is  analogous  to 
the  well-known  problem  of  reflection  and  transmission  of  transparent  sheets.  However, 
the  unusual  phase  relations  prevailing  in  the  case  of  metal-plate  media  make  certain  con¬ 
siderations  necessary  which  are  not  required  in  the  case  of  ordinary  dielectric  materials. 
For  this  reason,  a  theory  of  the  metal-plate  slabs  was  included  in  this  report.  It  is  devel¬ 
oped  along  the  lines  of  the  modern  theory  of  discontinuities  in  transmission  lines  introduced 
by  J.  Schwinger. 


The  problems  discussed  in  this  report  are  less  general  than  those  of  Carlson  and  Heins 
in  the  following  respect;  It  is  assumed  here  that  the  parallel  plates  are  so  stacked  that  the 
plane  containing  their  edges  is  perpendicular  to  the  plates,  while  Carlson  and  Heins  permit 
this  plane  to  make  any  angle  with  the  plates.  Perhaps  the  most  serious  limitation  of  this 
work  is  that  it  had  been  confined  to  a  plane  of  incidence  perpendicular  to  the  edges  of  the 
plates.  This  is  the  plane  of  the  maximum  constraint  from  the  point  of  view  of  wave  propa¬ 
gation  within  the  medium.  The  examination  of  other  planes  of  incidence  and  other  orien¬ 
tations  of  polarization  is  in  progress. 
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Terminology  and  Elementary  Properties 

A  parallel>plate  or  metal-plate  medium  is  a  part  of  space  containing  a  number  of 
parallel,  equidistant,  highly  conducting  sheets.  Such  a  structure  is  shown  in  Figure  1  with 
the  plates  perpendicular  to  the  plane  of  the  paper. 


It  is  assumed  that  the  plates  extend  to  infinity  forward,  back,  and  to  the  right,  that  the 
edges  of  the  plates  are  straight  and  perpendicular  to  the  paper,  that  the  plates  are  stacked 
Uniformly  as  shown  in  Figure  1.  Moreover,  to  avoid  the  complications  caused  by  edge  ef¬ 
fects,  it  is  assumed  that  there  are  many  plates  and  that  the  investigation  is  confined  to  the 
center  part  of  the  structure  so  that  one  plate  is  not  distinguished  physically  from  the  adja¬ 
cent  ones  by  virtue  of  its  position.  In  mathematical  parlance  this  assumption  is  expressed 
postulating  an  infinite  number  of  plates,  located  at  the  points  satisfying  the  conditions 
*  X  cotan  Q!,  X  =  ma,  m  =  0,  +1,  ±2,  ...  .  The  plane  through  the  edges  of  the  planes  is  called 
the  surface  of  the  medium. 

The  parameters  determining  the  electrical  properties  of  the  medium  are  the  plate  sep- 
^ation  a,  the  plate  thickness  t,  the  conductivity  a  and  the  shape  of  the  surface  of  the  medium 
<feter mined  by  the  angle  of  stacking  a.  The  propagation  of  electromagnetic  waves  within  the 
medium  is  determined  by  a  and  ct  ;  the  effect  of  the  latter  being  usually  negligible.  Surface 
phenomena,  such  as  reflection  of  an  incident  wave  will  depend  on  a  and  t  also. 

Metal-plate  media  are  generally  intended  for  use  in  a  well-defined  frequency  region 
mid  the  parameters  of  the  structure  are  so  chosen  that  t  is  very  small  compared  to  the 
^velength  X,  a  is  of  the  same  order  as  X  and  a  is  that  of  ordinary  metals  (5  x  10^  mhos 
meter),  in  first  approximation  one  always  assumes  t  =  0  and  a  =oo.  To  avoid  complications, 
study  is  restricted  to  the  case  a  =  90° . 
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Only  those  electromagnetic  waves  will  be  considered  whose  electric  vector  is  parallel 
to  the  metal  plates.  The  plate  separation  is  generally  so  chosen  that  at  the  given  frequeiJL 
only  one  mode  of  oscillation  is  propagating  in  the  medium.  This  condition  requires  ^ 


(1) 


It  is  well-known  that  there  is  no  propagation  within  the  medium  for  X  >  2a  and  that  for 
X  <  2a  the  ratio  of  the  phase  velocity  in  free  space  to  the  phase  velocity  in  the  medium  Is 


This  quantity  is  called  the  index  ctf  refraction  of  the  medium.  Let  k  be  the  propagation 
constant  in  free  space;  k  =  2ff/X,  and  let  x  =  nk  be  the  corresponding  quantity  in  the  medium. 
Equation  (2)  is  equivalent  to 


le  + 


(3) 


The  value  of  n  determines  the  propagation  within  the  plates  and  also  the  geometrical 
optics  of  the  rays  incident  on  the  surface  of  the  medium.  Rays  incident  in  a  plane  parallel 
to  that  of  the  plates  are  not  constrained,  while  rays  incident  in  other  planes  are  all  con¬ 
strained  by  the  plates.  The  classical  argument  pertaining  to  the  propagation  of  phase 
fronts  shows  that  unconstrained  rays  follow  Snell’s  law.  When  rays  are  incident  in  the 
plane  of  maximum  constraint,  the  (x,z)  plane,  the  Poynting  vector  in  the  medium  is  always 
directed  along  tBe  z-axls.  Still  one  might  speak  of  refraction  at  the  surface  in  terms  of 
equiphase  surfaces.  Consider  the  case  when  the  plates  are  not  staggered:  i.e.,  or  =  90P, 
and  let  the  radiation  be  incident  at  an  angle  e  as  shown  in  Figure  2.  The  points  Po,  Pi, 
P.1,  etc.  will  lie  on  a  discontinuous  equiphase  surface  which  is  approximated  by  the  plane 
shown  in  dotted  lines.  Let  the  Inclination  of  this  plane  be  Of.  Elementary  calculation 


gives 


sin  9 

tan  e. 
r 


=  n. 


(4) 


\ 

Fi6-  2  -  Refraction  in  the  plane  of  constraint 
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REFLECTION,  TRANSMISSION,  AND  THE  GRATING  EFFECT 

When  a  plane  wave  is  incident  on  the  surface  of  a  metal-plate  medium,  currents  are 
induced  on  the  plates.  In  addition  to  waves  set  up  in  the  medium  there  will  be  reflected  or 
scattered  waves  in  free  space.  Analogy  with  dielectric  materials  leads  one  to  expect  a 
transmitted  and  a  reflected  plane  wave.  However,  the  analogy  is  an  incomplete  one  since 
the  situation  at  the  surface  of  the  metal-plate  medium  differs  in  many  respects  from  that 
encountered  at  the  surface  of  dielectrics. 

First  of  all,  the  electric  field  at  the  surface  of  a  metal-plate  medium  is  not  uniform, 
but  a  periodic  function  of  position.  Therefore,  transmitted  and  reflected  waves  cannot  be 
uniform  plane  waves  as  they  are  in  the  case  of  dielectrics.  Only  at  a  distance  from  the 
surface  might  one  expect  approximately  uniform  plane  waves. 

There  is  never  more  than  one  reflected  wave  from  a  plane  dielectric  surface,  but  the 
edges  of  the  metal-plates  act  as  a  grating  and  for  a  single  plane  wave  incident  at  an  angle 
6,  there  will  be  reinforcement  of  scattered  waves  in  all  directions  0^  which  satisfy  the 
relation 

sin  fl  +  sin  =  m  ,  m  =  0,  *  1,  i  2,  .  .  .  .  (5) 


There  will  always  be  the  beam  6s  =  ~6  corresponding  to  m  =  0.  This  is  the  case  of  com¬ 
mon  (specular)  reflection  in  geometrical  optics.  There  might  be  other  beams,  provided 
that  toe  elements  of  the  grating  are  sufficiently  far  apart.  However,  when 

^  \ 

1  +  Isin0l<-  (6) 

then  (5)  cannot  be  satisfied  for  any  value  of  m  other  than  zero;  thus  there  is  only  one  re¬ 
flected  beam.  Since  a=  9(P  has  already  been  assumed,  one  may,  without  further  loss  of 
generality,  restrict  himself  to  positive  angles  of  incidence  and  omit  the  absolute  value 
signs  in  (6). 

The  electromagnetic  disturbance  caused  by  the  incident  plane  wave  at  the  surface  of 
the  medium  results  in  oscillations  in  the  space  between  the  plates.  The  latter  act  as 
Waveguides  of  uniform  width  and  infinite  height.  They  permit  the  passage  of  an  oscillation 
of  a  given  frequency  in  several  modes  propagating  with  different  velocities.  In  all  appli¬ 
cations  the  spacing  a  is  chosen  to  satisfy  the  inequality  (1),  so  that  in  the  contemplated 
frequency  range  only  the  fundamental  mode  is  propagated.  All  other  modes  excited  at  the 
o<lges  are  attenuated.  The  energy  associated  with  these  nonpropagating  modes  is  stored 
In  the  part  of  the  space  surrounding  the  edges.  The  presence  of  this  energy  is  responsi¬ 
ble  for  the  phase  shift  in  transmission  through  the  surface  of  the  metal-plate  medium.  If 
One  overlooks  the  field  present  in  these  higher  modes  he  is  led  tp  erroneous  boundary  con- 
<lltions  in  the  process  of  joining  the  incident,  reflected  and  propagated  waves  by  means  of 
^  continuity  argument.  TTiese  conditions  in  turn  lead  to  incorrect  reflection  and  trans- 
bJission  coefficients.® 

A  rigorous  theory  of  metal-plate  media  is  quite  difficult  and  has  not  yet  been  fully 
developed.  Carlson  and  Heins*  solved  the  problem  for  the  case  of  a  plane  wave  incident 
^the  plane  of  maximum  constraint  assuming  infinite  conductivity  and  zero  plate-thickness, 
^eir  mathematical  procedure  is  far  too  complicated  for  reproduction  here.  Only  the 
principal  results  will  be  stated  for  the  special  case  a  =  9(f, 


6 


NAVAL  RESEARCH  LABORATORY 


THE  RESULTS  OF  CARLSON  AND  HEINSf 

After  formulating  the  electromagnetic  field  problem  in  terms  of  an  integral  equattn 
of  the  Wiener.Hopf  type,  Carlson  and  Heins  proceed  to  solve  this  equation  by  the  use  of 
the  Fourier  transform.  The  crucial  step  in  the  solution  requires  the  splitting  of  a  comoi 
transcendental  function  in  two  parts;  one  of  which  has  no  zeros  or  poles  in  the  upper  h^^ 
plane,  the  other  is  similarly  behaved  in  the  lower  half-plane.  The  function  in  Question 
(forots90®)  “ 


l(w) 


_ _ sin^j^  -  w* _ 

-  w*|cos  a^^  -  -  cos  (ak  sin  9)j 


It  is  written  in  the  form 


(7) 


f{w)  = 


K.(w) 
K+(w)  * 


(8) 


where  K_(w)  has  no  poles  or  zeros  below  the  line  C  in  Figure  3  and  K^(w)  has  none  above. 


No  zeros  and  singularities  of  K^,  in  this  region 
I  (above  the  curve  C) 


fc 


Real  axis 


Fie.  3  -  The  plane  w 


■ffor  a  complete  understanding  of  the  next  two  sections,  "it  is  helpful  to  read  the  article 
of  Carlson  and  Heins,*  In  order  to  facilitate  a  comparative  reading,  the  convention  of 
using  a  time  dependence  in  the  form' of  e”^“*  has  been  retained  in  this  section.  Consequent¬ 
ly,  the  reflection  coefficient  is  written  in  the  form  of  r  •  |r|  .  To  eatvert  to  the 

standard  notation  used  in  other  parts  of  this  report  it  is  only  necessary  to  replace  i  by  -j. 


NAVAL  RESEARCH  LABORATORY 


7 


The  functions  are  properly  normalized  to  insure  desired  behavior  at  infinity.  The  actual 
form  of  K4.(w)  is 


.  V  _  1  w  ±.k  cosj  V(w)  Mg)  exD  (—  fin  2 
-  2  w  +  K  U(w)  I  TT 


(9) 


where 


with 


v(w)  =  1 

K-S) 

exp^ 

flaw 

;^2jrn 

ak  sin  9\ 

2itn  )• 

(10) 

w(w)  = 

exp 

/  Inw 
V  2vn 

ak  sin 

2Tn  / ' 

(11) 

\n  irn  / 

exp  • 

riawl 
l.vn  /' 

(12) 

4 

It 

L  ak  sin  5\* 
r  '  2jrn  / 

/ak  ' 
\2jrn, 

f. 

(13) 

II 

1 

/ak> 

r- 

(14) 

In  order  to  Insure  proper  behavior  of  U,  V  and  W  in  the  upper  half-plane,  the  sign  of 
and  e  ^  has  to  be  determined  as  follows:  When  -  0,  then  ^  0,  when  A^  <  0,  then 
0»  On  accoimt  of  the  restriction  expressed  by  (1),  is  always  positive  and  the  cor¬ 
rect  choice  is  E  Q  >  0.  The  numerator  of  (8)  then  takes  the  form: 


K.(w)  =-§ 


W  -  X 

w  -  k  cos  0 


Assuming  an  incident  field  of  the  form  Ex  =  0,  Ez  =  0,  Ey  =  =  exp  {ik(xsin  9  + 

*  cos  8)}  the  total  electric  field  is  represented  by  the  following  equation: 


«^(x,  z)  = 


I 


K^(k  cos  9)  I(w,  x)  e^^^  dw 


K^(w)  (w  -  k  cos  9)  sin  a  Vk^  -  w' 


T  > 


(16) 


^ere 


I(w,  x)=  [sin(x  -  n'a  -  a)  Vk* -  w*  _  gi^n.  sin  0  sin(x-n'a)  |4:®- w*}e^*'^^  ^  (17) 
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and  n'  is  the  largest  integer  not  exceeding  x/a.  This  is  the  special  case  (or  =  90®)  of  th* 
principal  result  of  Carlson  and  Heins  (p.  325).  The  path  of  integration  is  shown 

The  reflected  and  transmitted  fields  are  obtained  from  (16)  evaluated  for  negative 
positive  values  of  z  respectively.  One  can  show  that  the  behavior  of  the  integrand  in 
is  such  that  the  curve  C  can  be  closed  by  a  large  semicircle  in  the  upper  half-plane  for 
positive  values  of  z  and  by  a  similar  semicircle  in  the  lower  half -plane  for  negative  valu 
of  z.  The  contributions  from  these  circular  arcs  vanish  as  their  radii  tend  to  infinity, 
is  therefore  possible  to  evaluate  the  integral  in  (16)  by  means  of  the  theorem  of  residues 
the  poles  of  the  integrand  beii^  Wp  =  k  cos  6 ,  the  zeros  of  K+(w)  and  those  zeros  of  ’ 
sin  a  V^k*  -  w®  which  are  not  cancelled  by  the  numerator  of  the  integrand. 

In  the  case  of  the  transmitted  wave,  z  is  positive,  the  path  of  integration  is  closed  ia 
the  upper  half-plane,  therefore  only  poles  situated  above  C  need  to  be  considered.  Since 
K^.(w)  has  no  zeros  bi  this  region,  the  poles  are  at  Wn  =  k  cos  6  and  at  such  points  Wjjj  in 
the  upper  half-plane  which  satisfy  the  equation  a  Vk^  -  w^^  =  mx,  where  m  is  a  positive 
integer.  The  restriction  of  single  mode  propagation  within  the  medium,  expressed  by  the 
Inequality  (1),  has  the  consequence  that  all  Wm  with  m  >  1  are  pure  imaginary,  while 
Wi  =  (k®  -  77® /a*)  =  K.  For  m  >  1  it  is  then  proper  to  write  w^  =  iPm.  where  Pjjj  >  o. 

The  expression  obtained  by  evaluating  the  residue  at  Wm  will  contain  the  factor  e-PWz^ 
indicating  a  nonpropagating  wave.  Only  real  poles  generate  propagating  waves.  In  fact, 
the  residue  at  Wq  generates  the  one  at  w  =  x  generates  t  sin  (7rx/a),  where 


n'+  1  77  Kj.(wo)  1  +  e  ^  in'ak  sin  B  • 

Vic  K+(x)(Af  -  Wq)  ® 


(18) 


The  factor  si**  ®  takes  into  account  the  phase  delay  from  one  channel  of  the  medium 

to  the  next.  It  is  sufficient  to  concentrate  on  a  particular  pair  of  plates  and  to  examine 
only  the  field  between  them.  In  this  case  one  may  set  n'  =  0  and  simplify  the  expression 
of  the  transmission  coefficient  to: 


.  ...  J'T'  77  Kj.(k  cos  e) 

*  =  - va Kja)  — 


,  ika  sin  6 
1  +  e 


k  cos  B -  K 


(19) 


Assuming  that  all  deltas  defined  by  (13)  are  real,  one  obtains  after  a  lengthy  calculation:! 


-  2  ^  k  cos  B  _  2  ^  cos  B 
~  k  cos  B  +JC  ~  cos  0  + 1!  ’ 


(20) 


where  n  is  the  index  of  refraction  defined  in  (2). 

For  the  phase  of  the  transmission. coefficient  (19)  yields  readily: 

t'  a  k  sin  e  +  arg  K^(k  cos  B)  -  arg  K^(ic),  (21) 


f  Actually  Carlson  and  Htins  used  a  k  with  a  positive  imaginary  part,  integrated  along  the 
real  axis  and  then  Permitted  Im  k  to  approach  zero.  The  location  of  their  path  of  integra¬ 
tion  with  respect  to  the  singularities  is  identical  to  that  shown  in  figure 


t  See  Appendix,  Eqs.  (4 A)  and  (6AI. 
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provided  cos  0  >  n  which  is  usually  the  case.  The  first  term  in  this  expression  is  due  to 
the  location  of  the  origin.  It  disappears  when  the  origin  is  shifted  to  a  point  half-way  be¬ 
tween  the  edges,  such  as  Q  in  Figure  2. 

The  reflection  coefficient  is  obtained  by  evaluating  the  integral  in  (16)  for  large  nega¬ 
tive  values  of  z.  In  this  case  the  path  of  integration  is  closed  in  the  lower  half-plane. 

If  the  quantities  defined  by  (13)  are  all  real,  then,  apart  from  w  =  -k  cos  0,  K+(w) 
has  only  pure  imaginary  roots  in  the  lower  half-plane.  The  latter  generate  attenuated  waves 
as  was  shown  before.  The  roots  of  sin  a  Vk®  -  vT  in  the  lower  half-plane  are  cancelled  by 
the  poles  of  K+(w),  therefore  there  is  one  propagating  reflected  wave.  It  is  of  the  form: 


(22) 

(23) 


where  is  dK+/dw.  Equation  (23)  gives: 

p>  =  arg  K^(k  cos  0)  -  arg  K'^(-k  cos  0)  +  v.  (24) 

When  the  deltas  in  the  functions  V  and  W  are  all  real,  I  rl  can  be  computed  by  elemen¬ 
tary  means.  In  this  case:t 


COS  0)  ^  2k  cos  0  I  cos  0  t  K  I 
^(-k  COS0)  cos  y  I  ^  0  4.ii  I’ 


and  hence  from  k=  nk  and  equation  (23)  follows 


**''cos0+n 

The  absolute  value  sign  has  been  omitted  on  the  right  because  the  restriction  on  the 
•deltas  implies  cos  0  >  n  .  This  will  be  shown  Immediately.  First,  it  is  proper  to  examine 
meaning  of  the  restrictions  ^  0  for  all  values  of  n.  For  negative  values  of  n,  A^  is 
positive  for  every  angle  of  incidence  0  (in  the  first  quadrant).  This  is  a  consequence  of 
coi^tion  (1),  imposed  to  insure  single  mode  propagation  between  the  plates.  The  condition 
^n  =  0  for  positive  n  is,  however,  more  restrictive.  Since  A^  <  An+i .  needs  to  be 

examined.  The  latter  is  not  negative  if,  and  only  if. 

Or 

X  2v  >  '  ,  -  (28) 

_ _ 

Iw*  doctors  resulting  from  the  infinite  products  cancel  out,  since  for  any  real  number  u  , 

~  uhen  all  deltas  are  real  and  |  |  =  \U(wa)\  be- 

*«  the  ePsilons  are  always  real  on  account  of  condition  (i).  The  evaluation  of  in  the 
’**ral  case  is  discussed  in  the  Appendix, 


I 

r  - 
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corresponds  to  an  angle  of  incidence  0  for  which  A*  <  0,  but  A|  >  0.  According  to- (13) 

Clearly  -A?  <(ak/2jT)*,  therefore  whenever  A*  <  0  there  is  an  angle  6'  such  that 

Ai  =  -i  ^  cos  6'.  (33) 

From  (32)  then  follows 

±sin0'=^  -sinfl.  (34) 

It  is  proper  to  restrict  oneself  to  the  use  of  the  +  sign  in  (34)  since  (33)  does  not  specify 
the  sign  of  6'.  Then 

sin  6  +  sin  ^ 


The  function  K+(w)  now  has  an  additional  real  root,  since  the  factor  Ai  -iaw/2jr  van¬ 
ishes  when  w  =  -k  cos  0' .  When  the  presence  of  this  real  pole  is  taken  into  account  in  the 
evaluation  of  the  Integral  in  (16),  a  new  propagating  wave  is  found  in  the  region  z  <  0.  In 
fact,  from  (17)  and  (35)  follows  I(-k  cos  6',  x)  =  -sin(ak  sin  0')  exp(-ikx  sin  0')  and  the 
contribution  of  the  pole  is  exp{4k(x  sin  6'  +  z  cos  6'},  where 

(1)  K^(k  cos  0) _ 1 _  ^3gj 

^  ~~  K'^(-k  cos  0')  k(cos  0  +  cos  0')  * 

This  is  the  first  order  diffracted  beam;  it  propagates  at  an  angle  0'  with  the  surface  normal. 

The  quantities  necessary  for  the  calculation  of  |tl,  jrj  and  |r*’  |  are  evaluated  in  the 
Appendix  for  the  case  when  A?  ^  0,  but  >  0.  The  following  e:q)ressions  result 


i'ti' 

■f 


t 


!. 


f 


2 cos  0  i/cos  0'  -  n  cos  0’  +  cos  0  ,  (37) 

cos  e  +  n  f  cos  0'  +  n  cos  0'  -  cos  0 


cos  8  -  n  cos  0  +  cos 
cos  0  +  n  cos  0  -  cos  0' 


(38) 


2  cos  0 

|cos  0'  -  cos  01 


]/ 


cos  0  n 
cos  0  +  n 


n  -  cos  0'  , 
n+  cos  0' 


(39) 
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These  equations  are  valid  for  those  angles  of  Incidence  6  which  satisfy  the  inequalitig 

They  can  be  put  in  a  symmetrical  form  by  the  introduction  of  the  angle  Oq  =  arc  cos  n. 

It  is  Interesting  to  examine  these  coefficients  for  the  angle  of  incidence  6  =  Qq.  Pro^, 
cos  ©o  =  n  follows  sin  Oq  =  x/2a.  Then  (35)  takes  the  form: 

sin  e  +  sin  0'  =  2  sin 

0  (41) 

Therefore,  if  two  of  the  angles  0,  O’,  Oq  are  equal,  then  all  the  three  are  equal.  The  ex¬ 
pressions  (37),  (38),  and  (39)  are  of  the  indeterminate  form  0/0  for  0  =  Oq. 

On  evaluating  lim  cos  0  -  cos  ^  according  to  L’Hospital's  rule  and  noting  that 
A  ^  COS  “  “  cos  t/ 

C^t70 

dO’/dO  =  -1  when  0  =  O',  one  obtains  |r|  =  1/2  for  0  =  Oq.  By  comparing  (23)  and  (36)  it  is 
clear  that  |r  |  and  |r^  |  must  be  equal  for  0  =  O’,  therefore  [r^*  |  =  1/2.  A  brief  calculation 
yields  jt|  =  1. 


ENERGY  RELATIONS 

Since  the  experiments  are  generally  so  designed  that  energy  flow  and  not  amplitude  is 
detected,  it  is  of  practical  value  to  derive  expressions  for  the  former.  The  quantity 
Pr  =  Ir  |*is  a  correct  measure  of  the  relative  power  in  the  beam  reflected  according  to  the 
law  of  geometrical  optics.  The  quantity  1 1 1  is  the  maximum  of  the  electric  vector  between 
a  pair  of  plates,  the  space  variation  being  sinusoidal.  Thus  the  average  rate  of  energy 
flow  in  the  direction  of  the  z-axis  is  Inltf*.  The  change  in  cross  section  of  the  beam  at  the 
surface  Introduces  a  factor  sec  0,  therefore  the  relative  transmitted  power  is: 


=  »  Hi  1 1*  sec  0  = 


4  n  cos  0 
{n  + cos  of 


for  the  range  specified  by  (28)  and 

_  4  R  cos  0  cos  0’  -n  cos  O'  +  cos  0 
^t  “  (n  +  cos  o  f  cos  0'  +  n  cos  0’  -  cos  0 


(43) 


for  the  range  specified  by  (40). 

In  general,  there  is  a  change  in  cross  section  in  the  case  of  the  diffracted  beam. 
Therefore 


cos  0’ 
cos  0 


(44) 


when  0  satisfies  (40)  and  Pr‘“  =  0  when  0  satisfies  (28). 
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The  quantities  Pj.,  Pr^^^and  p^  describe  the  division  of  power  incident  on  a  surface 
element  of  an  air -metal-plate  medium  interface  of  infinite  extent.  The  conservation  of 
energy  requires  pj.  +  pj.'v  +  p^  =  1.  Figure  5  shows  the  variation  of  power  division  with 
angle  of  incidence  for  2a/x  =  1.10,  1.20,  1.30  and  1.40. 

When  6  =  6’  =  6 Q,  the  Pj.  =  p^^^)  =  1/4  and  pt  =  1/2;  i.e.,  one-half  of  the  energy  is  trans¬ 
mitted  into  the  metal-plate  memum,  the  other  hall  is  equally  divided  between  two  beams 
located  symmetrically  about  the  surface  normal. 

L.  J.  Chut  previously  deduced  this  division  of  power  for  0  =  0o  by  elementary  con¬ 
siderations.  He  wrote  the  incident  wave  in  the  form 


Each  term  on  the  right  of  (45)  corresponds  to  a  wave  carrying  one-half  of  the  energy.  The 
second  term  satisfies  the  boundary  conditions  within  the  metal-plate  medium;  it  is  the 
transmitted  wave.  The  first  term  represents  a  wave  antisymmetric  about  the  plane  x  =  a/2; 
it  is  responsible  for  the  storage  of  energy  at  the  interface  and  for  the  presence  of  reflected 
waves.  By  elementary  symmetry  considerations  one  can  deduce  that  the  reflected  waves 
make  equal  angles  with  toe  surface  normal  and  carry  equal  amounts  of  power. 

It  is  of  some  interest  to  note  that  if  toe  index  of  refraction  of  a  metal-plate  structure 
is  chosen  large  for  toe  purpose  of  reducing  reflection  at  normal  incidence,  the  advantage 
so  gained  is  quickly  lost  when  toe  same  structure  is  illuminated  at  an  angle  which  per¬ 
mits  toe  appearance  of  a  diffracted  beam.  This  is  born  out  by  Figure  5. 


media  bounded  by  two  parallel  surfaces 

The  theory  described  in  the  preceding  section  gives  toe  reflection  and  the  transmission 
coefficients  of  the  air-metal>plate  interface.  In  any  experimental  test  in  which  the  energy 
propagating  between  the  plates  is  not  completely  absorbed,  it  is  necessary  to  take  into 
account  the  multiple  reflections  within  the  parallel-plate  structure.  A  metal-plate  medium 
of  finite  uniform  depth  bounded  by  plane  surfaces  will  be  called  a  metal-plate  slab.  The 
reflection  from  such  a  slab  depends  on  the  reflection  coefficients  of  the  front  arid  of  the 
back  sxirfaces  and  on  the  distance  between  these  surfaces. 

From  the  point  of  view  of  theory,  the  basic  quantities  are  the  coefficients  describing 
the  air-metal-plate  interface.  The  experimenter,  however,  always  operates  with  slabs  or 
other  shapes  of  fifaite  depth.  It  is  possible  to  reduce  greatly  toe  effect  of  multiple  reflec¬ 
tions  by  constructing  deep  slabs  containing  absorbing  wedges,  but  some  uncertainty  with 
regard  to  toe  phase  and  magnitude  of  the  back  surface  reflections  still  remains.  It  is 
therefore  desirable  to  deduce  expressions  describing  the  behavior  of  slabs,  or  metal-plate 
®edia  bounded  by  two  parallel  surfaces. 

The  corresponding  problem  has  been  solved  in  connection  with  dielectric  sheets  and 
lossless  transmission  lines.  Let  n  be  the  reflection  coefficient  and  L  the  transmission 

t  Oral  communication;  see  Acknouledgments. 
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iBBBBBBBBBBBBBBBBBBBBBK'BBBBBBl 

iBBBBBBBBBBBBBBBBBBBBBiaKBBBBBi 

laaBBBBBBBDXa^IlZESSrElBBBRBBI^BBBBl 

iBBBBBBBBBUQBaBBBBBBBaBBB^BBBl 

iaBBBBBBBBBBBaBBBBBBaBBBBaiSBBi 

iBBBBBBBBBBBBBnBBBB^BBBBBBBral 

iBBBBBBBBBBBBBBKilKaBfiSaBBBB^I 

iBBBBBBBBBBBBBiSaBaaBaBBBliSiBBBi 

iBBBCriTiTnr^BBBKBBBaBBBBBBBBBBBIil 
|BBBBBBB!IBBB'JBB^BBBBBBBBBBBBB| 
- IBBBBBB9r;BBl!;<BBBBBBBaBBBBBB| 

iBBBBaaBiisiiisir 


1.30,  n  -  .639 


— •  1.40,  n 


ig.  5  -  Power  distribution  as  a  function  of  angle  of  incidence 
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coefficient  of  the  interface  separating  media  1  and  2  when  the  radiation  is  incident  from 
medium  1  (usually  air)  and  let  ft  and  t2  have  similar  meanings  for  a  radiation  incident 
from  medium  2.  These  quantities  will  in  general  depend  on  the  angle  of  incidence  6.  To 
distinguish  coefficients  pertaining  to  the  slab  from  those  pertaining  to  single  surfaces, 
the  former  will  be  denoted  by  capital  letters.  One  can  show  that  the  reflection  coefficient 
of  a  slab  of  medium  2  surrounded  on  both  sides  by  medium  1  is  given  by  the  equation! 


where  6^  is  the  phase  delay  between  the  two  intersections  of  a  normal  of  the  slab  with  the 
surfaces.  In  the  case  of  a  lossless  dielectric  slab  of  thickness  d  and  of  refractive  index 
n  one  has  62  =  dk2  cos  62,  where  k2  =  nki  =  2tT\/\  and  Q2  is  the  angle  of  refraction. 

Similarly,  for  the  transmission  coefficient  one  finds  in  the  literature 


A  few  words  of  caution  about  the  phases  of  these  coefficients  are  in  order.  The  coef¬ 
ficient  R  refers  to  the  front  surface  of  the  slab.  If  the  plane  of  reference  is  chosen  else¬ 
where,  say  in  the  plane  of  symmetry  of  the  slab,  the  phase  of  R  must  be  altered.  Equation 
(47)  gives  the  phase  of  the  transmitted  wave  at  (x,  y,  d  +  0)  with  respect  to  the  phase  of  the 
incident  wave  at  (x,  y,  0),  and  not  the  change  of  phase  produced  by  the  introduction  of  the 
slab  into  the  path  of  the  radiation.  If  the  latter  is  desired,  one  has  to  multiply  Tiby 
exp  jCi,  where  61  is  the  phase  delay  through  a  slab  of  thickness  d  made  of  medium  1.  When 
medium  1  is  air,  which  is  usually  the  case,  61  =  2ird  cos  0/X. 

Let  the  incident  wave  be  of  the  form  exp  {-jk(z  cos  6  +  sin  0)}’  Jand  let  the  slab  be 
bounded  by  the  planes  z  =  0  and  z  =  d.  The  reflected  wave  will  then  have  the  form 
R  expj[-jk(-z  cos  0  +  x  sin  0)J’and  the  transmitted  wave  in  the  region  z  >  d  will  have  the 
form  TTexp  {^jk(z  cos  0  +  x  sin  0)J-,  where  T  =  Ti  ei^.  For  slabs  bounded  on  both  sides 
by  the  same  medium,  the  coefficient  T  seems  more  fundamental  than  Tj,  although  the  latter 
is  more  commonly  found  in  the  literature.  T  does  not  depend  on  a  plane  of  reference;  it 
compares  the  wave  in  the  presence  of  the  slab  with  the  wave  in  its  absence. 

The  principal  difference  between  dielectric  slabs  and  metal-plate  slabs  is  that  in  the 
case  of  the  latter  the  single  surface  coefficients  rj,  Tz,  t^  and  t2  are  not  real.  The  simple 
equations  connecting  these  quantities  must  be  replaced  by  new  ones.  These  will  now  be 
derived  for  the  case  when  the  angle  of  incidence  is  such  that  there  is  only  one  reflected 
beam.  The  discussion  can  be  made  general  so  as  to  include  lossless  dielectric  slabs,  and 
discontinuities  in  transmission  lines  as  well. 


general  discontinuity  relations 


Let  the  plane  z  =  0  be  a  locus  of  discontinuity  for  the  propagation  of  electromagnetic 
waves.  This  plane  might  be  the  boimdary  of  two  dielectric  media,  the  surface  of  a 

tRef.  1,  p.  355. 

tihe  time  dependence  is  exp  jojt. 
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metal-plate  medium,  a  place  where  the  dimensions  of  a  transmission  line  undergo  a  sudden 
change,  etc.  In  general,  the  propagation  constants  on  the  two  sides  of  the  discontinuity  win 
be  different. 

To  simplify  matters  it  is  practical  at  this  point  to  exclude  multiple  propagation  and 
multiple  reflection.  While  in  the  neighborhood  of  the  discontinuity  energy  might  be  stored 
in  the  form  of  “higher"  modes, it  is  assumed  that  far  away  from  the  discontinuity  on  either 
side  there  is  only  one  wave  propagating  toward  and  one  away  from  the  discontinuity.  The 
case  of  plane  waves  arriving  at  the  plane  z  =  0  at  oblique  angles  of  incidence  is  included 
in  this  discussion.  To  each  angle  of  incidence  6^  in  the  medium  to  the  left  of  the  discon¬ 
tinuity  the  law  of  refraction  orders  an  angle  flj  in  the  medium  to  the  right.  The  four  waves 
to  be  considered  together  have  the  directions  of  propagation  flj,  ir -  ^  ir  -  4  .  These 
are  shown  in  Figure  6. 


It  is  appropriate  then  to  observe  the  discontinuity  from  far  away  on  both  sides,  so  far 
that  one  has  to  deal  with  a  pair  of  waves  on  each  side.  Formally  such  a  system  is  identical 
with  a  two-terminal  pair  network.  The  plane  (or  propagating)  waves  observed  far  from  the 
discontinuity  can  then  be  extended  mathematically  to  the  plame  z  =  0,  and  one  might  speak 
of  the  amplitude  of  the  plane  wave  at  z  =  0  even  though  a  very  complex  field  exists  in  the 
latter  neighborhood.  When  for  large  negative  values  of  z  one  has  a  wave  of  the  asymptotic 
form  bj  exp  ■fjkz},  the  quantity  b^  will  be  regarded  as  the  amplitude  of  this  wave  at  z  =  0. 
This  is  the  customary  procedure  employed  in  discussion  of  discontinuities  in  waveguides.* 

Let  a^  and  a^  be  complex  numbers  proportional  to  the  electric  amplitudes  at  z  =  0  of  the 
waves  traveling  toward  the  discontinuity  and  let  b^  and  b^  have  a  similar  meaning  for  waves 
traveling  away  from  the  discontinuity.  The  subscripts  are  indicative  of  the  media  1  (left) 
and  2  (right).  Let  a  be  a  vector  with  components  aj  and  aj  and  let  b  have  a  similar  meaning. 
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On  account  of  the  linear  character  of  Maxwell’s  equations  and  the  boundary  conditions  the 


vector  b  depends  on  a  linearly: 


b  aSa. 


(48) 


The  scatterii^  matrix  8  connects  the  incoming  and  the  outgoing  waves.  It  is  easy  to 
see  that 


It  has  been  shown  that  when  the  vectors  a  and  b  are  properly  normalized,  the  matrix 
$  is  symmetrical.t  - 

When  no  dissipation  of  energy  takes  place  at  the  discontinuity  the  conjugate  matrix  8* 
Is  the  Inverse  of  S.J:  In  symbols: 

SS*  =  I,  (50) 

where  I  is  the  identity  matrix.! 


The  matrix  equation  (50)  is  equivalent  to  the  following  four  scalar  equations: 


\ 

1'/ 

(51) 

Irj-.t.t.M, 

(52) 

rt*+tr**0, 
a  2  a  ’ 

(53) 

From  (51)  and  (52)  follows: 

t  r*  +  r  t*  =0. 

.11  a  1 

154) 

-Faf  =tiV-ti*tg=2Imt^V. 

(55) 

G.  Ref.  4,  i>.  148.  A  matrix  S  called  symmetrical  uHen  it  is  identical  with  the 
^’’^sfiosed  matrix  S  obtained  by  interchanging  rows  0/  S  with  its  columns.  Proper 
*^'^ali zation  of  the  vectors  and  b  means  that  the  factors  of  proPortioi^ality  ietween 
t^»r  components  and  the  electrical  quantities  are  so  chosen  that  i  lojl  ,  i  IojI  ^ 

and  jf  jisl*  represent  the  mean  energy  flows  in  the  four  waves  involved.  The 
^“lues  of  rj  and  r^  are  indSpendent  of  the  normalization  and  so  is  the  product  tj  tg 
Since  (46)  and  (4T)  contain  tj,  and  tg  only  in  the  form  of  their  product,  so  far  as  slabs 
Concerned,  the  normalization  is  quite  immaterial , 

tthe  symbol  *  indicates  the  complex  conjugate  of  the  quantity  it  follows. 

^9^*  Ref.  4,  p.  14Q.  A  direct  proof  of  (so)  can  be  had  by  reversing  the  direction  of 
and  thus  interchanging  incoming  and  outgoing  waves.  This  change  also  replaces 
*^ery  element  of  S  by  its  canplex  conjugate,  since  time  reversal  is  equivalent  to  the 
^*fiacement  of  j  by  -j.  One  obtains  in  this  way  a  which  togethetr  with  (48)  suffices 

”  establish  (so). 
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The  equality  of  the  real  number  on  the  left  of  (55)  with  a  pure  Imaginary  on  the  right  in, 
plies  that  they  are  both  zero,  therefore 


\^i\  =  \h\ 


(5«; 


and  titz*  is  real.  In  fact,  from  (52)  follows:  hta*  =  1  -  -  0  because  |r2 1  cannot  exc 

one.  Therefore  arg  ^  =  arg  Let  p  =  1  ri  1  and  t  =  yiiia’*  =  Vl  -p*.  One  may  write: 


r 

1 


.-jp" 


jr' 


•^  =  pe  ,t^  =qTe  "  ,  t^  =  q 


-1  -j 

—  n  * 


re 


(57) 


where  the  positive  factor  q  depends  only  on  the  normalizatlonf  of  the  amplitudes.  From 
(53)  then  follows: 

exp  K-p'  +t')=  -  exp  j(p"  -  T')  =  exp  j(p"  -  t'+  x).  /g. 


Hence: 

p' +  p"  -  2T'=±7r.  (59J 

This  last  relation  is  of  fundamental  importance.  In  the  case  of  a  dielectric  sheet  in  air 
p '  =  jr,  p  "=  0  and  t'=  0.  In  the  case  of  metal-plate  media  p ',  p"  andr'  depend  on  the 
frequency. 

It  has  been  assumed  so  far  that  the  origin  is  located  at  the  discontinuity.  This  condi. 
tion  will  now  be  removed  and  the  concept  of  discontinuity  generalized  in  preparation  for 
the  discussion  of  slabs. 

The  position  of  the  origin  is  physically  immaterial,  a  shift  of  the  origin  leaves  the 
magnitudes  of  the  reflection  and  transmission  coefficients  unchanged.  If  a  new  coordinate 
system  is  introduced  with  its  origin  O'  at  (0,  0,  1),  then  ai  exp(-jki  z  cos0 )  and  aa  exp 
(jkj  z  cos  02 )  will  be  replaced  by  ai'  exp  (-jkiz'  cos  6)  and  as'  exp(jk2z'^os  02 )  respectively, 
where,  on  account  of  z  =  z'  +  1 

_  J(A  ,  jA 

a'  =  aj  e  ■’^x,  a^'  =  agO'^a 

with  A 1  =  ki  1  cos  0,  As  =  k2l  cos  02 .  t  Similarly  bi  and  bs  are  replaced  by  bi'  =  bi  e^^i  and 
bs' =  bse'jAa  respectively.  ' 

The  relations  between  the  amplitudes  in  the  two  systems  of  reference  can  be  expressed 
in  a  compact  form  by  the  use  of  the  matrix 


t  See  footnote  on  p.  I’j.  > 

t  When  this  theory  is  applied  to  discontinuities  in  transmission  lines  0  ^  «  o  and  ki  and 

k2are  the  pertinent  phase  constants  of  propagation . 
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In  fact 

&'  =  F'a,  b'  =  Pb. 

(61) 

Then  (48)  becomes 

b'  =  PSPa', 

(62) 

showing  that  the  shift  of  origin  changed  the  scattering  matrix  S  into  PSP.  Elementary 
matrix  calculation  results  in: 


PSP 


/riexp2j(^i  tjexpj((#^ 

\  ti  exp  j(<^^  r^  exp  -  2j0j  /  * 


(63) 


In  order  that  the  preceding  results  be  valid,  it  is  not  necessary  that  there  should  be  a 
single  plane  or  point  of  discontinuity.  It  is  sufficient,  that  in  some  part  of  the  space,  say, 
for  large  negative  z,  there  should  be  two  plane  waves  progressing  in  opposite  directions 
(or  in  directions  malting  aisles  6i  and  n  -  9i  with  the  z-axis),  and  that  in  another  part  of 
the  space  there  should  be  a  similar  pair  of  waves.  The  amplitudes  of  these  waves  should 
be  linearly  related  and  there  must  be  no  dissipation  of  energy  so  that  every  process  is 
time  reversible.  It  is  clear  then  that  lossless  discontinuities  located  in  lossless  media 
can  be  lumped  together  into  a  single  lossless  discontinuity  in  the  same  way  as  a  series  of 
four  terminal  networks  can  be  lumped  together  into  a  single  four  terminal  network.  The 
single  discontinuity  representing  the  entire  system  of  discontinuities  can  then  be  shifted  to 
any  convenient  point.  The  matrices  associated  with  a  location  of  this  composite  discontin¬ 
uity  are  the  elements  of  the  matrix  family  PSP,  where  S  is  the  matrix  associated  with  a 
fixed  location  of  the  discontinuity  with  respect  to  the  coordinate  system. 

The  quantities  <t>i  and  entering  into  the  definition  of  P  have  to  be  determined  with 
some  care  since  it  is  not  at  all  clear  (or  relevant)  what  happens  in  the  neighborhood  of 
the  origin  where  the  discontinuities  might  be  located.  A  sMft  of  the  origin  must  be  thought 
of  as  a  shift  of  the  reference  system  of  the  observer,  therefore  in  the  equations  ifn  =  ki  Icos  9 
and  ^2  =  k2 1  cos  92  the  appropriate  values  of  k  and  9  are  those  which  prevail  at  the  site  of  ob¬ 
servation;  i.e.,  far  removed  from  the  discontinuity.  When,  for  example,  the  composite  dis¬ 
continuity  consists  of  a  dielectric,  or  metal-plate  slab  surrounded  by  air,  then  ^  =  0%  since 
the  observation  occurs  in  air,  k2  =  ki,  and  the  emergent  beam  is  parallel  to  the  incident  one. 
This  situation  prevails  for  a  large  class  of  discontinuities  which  deserve  closer  examination. 

When  a  discontinuity  possesses  a  high  degree  of  symmetry  this  fact  is  reflected  in  the 
natrix  S  provided  that  Uie  system  of  reference  is  properly  chosen.  In  the  case  of  a  slab 
(dielectric,  or  metal-plate)  or  a  symmetric  iris  in  a  waveguide,  the  plane  z  =  0  may  be 
made  to  coincide  with  the  plane  of  symmetry  of  the  structure.  A  reversal  of  the  positive 
<llrection  of  z  will  then  leave  the  matrix  S  unchanged.  This  reversal  is  equivalent  to  an  inter¬ 
change  of  the  indices  1  and  2.  Therefore  ri  =  r2,  ti  =  t2.  From  (59)  then  follows: 

|P'-T'|=|.  (64) 

This  result  is  formalized  as  a  theorem. 

phases  of  the  reflection  and  the  transmission  coefficients  of  a  symmetric  lossless 
discontinuity  differ  by  90^  provided  that  the  point  of  reference  is  the  center  of  symmetry. 


In  experimental  work  it  is  often  convenient  to  refer  the  phases  of  the  waves  involved 
Jo  one  of  the  physical  surfaces  of  the  discontinuity  such  as  the  first  air-dielectric  interface 
the  case  of  a  dielectric  sheet.  Let  Sq  be  the  matrix  associated  with  the  symmetric  dis¬ 
continuity  with  the  origin  located  in  the  plane  of  symmetry  and  let  d  be  the  over-all  thick- 
hoss  of  the  structure.  When  the  origin  is  shifted  to  the  front  surface,  1  =  -d/2,  the  matrix 
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So  is  replaced  by  S  =  PSqP,  with  (h  =4>  2  =  -^kdcos  6,  where  k  is  the  propagation  constan* 
of  the  surrounding  medium.  From  (63)  follows: 


therefore r  =  rQ  and  p'  =  Pg'  -  20i  =  p©  +  dkcos  fti.  It  is  clear  that  the  quantities  p,  t  andT' 
have  a  physical  significance  independent  of  the  location  of  the  reference  plane.  This  fact 
manifests  Itself  mathematically  by  the  invariance  of  these  quantities  under  the  transfor¬ 
mation  PSP  with  =  ^2. 

A  lossless  symmetric  discontinuity  is  electrically  characterized  by  two  quantities 
say,  pandT',  since  t  =  \f\  -  p^  and  p'has  to  satisfy  (64). t  The  meaning  of  two  electrical 
parameters  is  easily  understood.  The  physical  parameters  available  for  adjustment  are 
the  series  Impedance  and  the  shunt  reactance.  In  the  familiar  case  of  a  dielectric  slab 
there  is  no  energy  storage,  therefore  the  shunt  reactance  in  the  equivalent  circuit  is  zero. 
This  is,  however,  not  the  case  for  metal-plate  slabs. 


FORMULAS  FOR  METAL-PLATE  SLABS 

A  slab  of  metal-plate  medium  surrounded  by  air  is  a  symmetric,  approximately  loss¬ 
less  discontinuity  and  the  results  of  the  last  few  paragraphs  apply.  In  order  to  test  the 
predictions  of  the  Carlson  and  Heins  theory  it  is  sufficient  to  determine  experimentally  the 
magnitude  of  the  reflection  coefficient  and  the  phase  of  either  the  transmission  or  the  re¬ 
flection  coefficient. 

With  the  aid  of  equations  (57)  and  (59)  one  can  derive  from  (46)  by  elementary  manipu¬ 
lations 


^(1  -  p^*+  4p*sin*  \l/ 


where  i(/=  p"  +  62.  The  absence  of  loss  immediately  justifies  the  equation 

(67) 

(68) 

it  follows  that 

T'  =  2t'-  +  6^  +  arg  -  p*  exp  [-2j(p''  + 

Actually  there  is  an  uncertainty  in  sign  resulting  in  an  uncertainty  of  rr  in  the  phase 
of  the  reflection  coefficient. 
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But  1 

and  therefore 

-  p®exp  [-  2jiJ/]  =  - 

^(1  -  p*)  cosi//  +  j(l  +  p*)  sin^  exp  [-  lip] 

where 

arg<{ 

1  -  p*  exp  2j  =  i//'  -  }j/, 

(70) 

=  arc  tan  ^  tan  i//\ 

(71) 

and  the  multiple  valued  function  arc  tan  x  is  so  chosen  that  i^and  in  the  same  quad¬ 
rant.  From  the  last  three  equations  and  equation  (59)  it  follows  that 


T'  =  p'  + ^'+ ir  -  fij.  (72) 

The  phase  of  R,  if  required,  can  now  be  computed  with  the  aid  of  (64)  which  applies  in  this 
case  to  the  quantities  denoted  by  capital  letters. 


AMPLITUDE  COMPUTATIONS 

The  numerical  computation  of  Ir  |,  Ul  and  |i^') I  presents  no  difficulties.  First  the 
limiting  angle  6-^  is  foundfrom  Figure  4,  For  6  ^  6l>  I*"  I  found  from  (26),  |tl  from  (20) 
and  =  0.  For  9  >  ^l,  equations  (37),  (38)  and  (39)  apply.  Some  care  has  to  be  exer¬ 
cised  in  the  neighborhood  of  the  point  0  =  0  o»  because  all  fractions  involved  are  of  the 
form  0/0.  This  indeterminacy  is  easily  removed  by  cancelling  factors  vanishing  for 
9  =  Aq.  With  the  aid  of  (41)  one  may  transform  (38)  to  the  form: 


Ir  1  = 


1  (COS9'.COg,tf 
4  cos  ^0+  cos  9 


9  +  9o 
2 


cosec  do 


(73) 


The  same  procedure  is  applicable  to  |r(‘)|  and  |t  |.  Figure  7  shows  the  values  of  p-=  |r| 
plotted  against  9  for  selected  fixed  values  of  x  =  2aA .  Table  1  gives  the  same  information 
in  numerical  form  and  in  greater  detail.  While  the  shape  of  the  curve  p[9)  resembles  that 
of  the  reflection  coefficient  of  dielectric  media  for  parallel  polarization,  the  curves  here 
plotted  do  not  reach  the  axis  p  =  0. 


i; 


PHASE  COMPUTATIONS 

The  numerical  computation  of  the  phases  of  the  reflection  and  transmission  coefficients 
is  a  lengthy,  difficult  task,  even  in  the  simplest  case  when  a  diffracted  beam  is  not  present. 
On  accoimt  of  limitations  in  manpower  assigned  to  this  project  only  those  coefficients  were 
computed  which  have  a  direct  connection  with  the  experimental  program  undertaken.  The 
results  appear  at  the  end  of  this  report. 

Since  the  quantities  p' ,  p"  and  r'  satisfy  (59),  it  is  sufficient  to  compute  two  of  these. 
Only  p"  enters  into  the  expression  for  1r  t  the  magnitude  of  the  reflection  from  a  slab  of 
metal -plate  material.  This  is  the  quantity  most  easily  accessible  to  measurement  and  for 
this  reason  the  major  part  of  the  computing  effort  has  been  devoted  to  the  computation  of 
p"  as  a  function  of  x  =  2a/x  and  9 .  For  normal  incidence  p'  was  also  computed  as  a  func¬ 
tion  of  X.  This  quantity  will  be  considered  first.  Next  p"  will  be  calculate  for  normal 
incidence,  finally  p"  for  oblique  Incidence.  In  this  last  case  the  angle  of  incidence  will  be 
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0  in  dn^recs 

Fig.  7  -  Magnitude  of  the  reflection  coefficient  as  a 
function  of  angle  of  incidence  (specular  reflection) 


SO  restricted  that  no  diffracted  beam  occurs.  The  value  of  p'  Is  given  by  equation  (24). 

In  view  of  equation  (9)  this  can  bp  given  the  following  form: 

argK^'  (-kcos0)  =  argV(-kcos6)  +  argW(-kcos0)-  argU(-kcose)-arg(-kcose  +  ic) 

-  ak  cos  0  (In  2  -  1)/^.  (74) 

When  there  is  only  one  reflected  beam  the  fimctionsV(-x),W(-x)  and U(-x)  are  complex 
conjugatesofthefunctlons  V(x),W(x)andU(x)respectively,  moreover  -k  cos  9  +  K  is  negative, 
therefore 

arg  K^'  (-k  cos  6)  =  -arg  (k  cos  6)  +  v.  (75) 

When  there  are  two  reflected  beams  (75)  is  again  valid,  but  one  has  to  verify  it  sepa¬ 
rately  for  6>  6o  and  for  0  <  ft,.  Consequently, 

p'  =  2  arg  (k  cos  9).  (76) 

It  is  now  necessary  to  calculate  arg  K^.(k  cos  0)  for  the  case  0  =  0.  In  this  case  from 
(13)  and  (14)  follows:  _ _ 
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therefore  V(w)  =  W(w).  Also  |  An  +  i(ak/2jrn)  |*  =  1. 

It  is  practical  to  introduce  the  variable  x  defined  as  follows: 

_  ka  ^  (78) 

TT  ~  X  • 

The  values  of  x  are  limited  to  the  interval  1  to  2.  For  practical  reasons  the  computations 
are  confined  to  the  range  1.01  to  1.35.  From  (9)  then  follows: 

arg  (k)  =  2arg  V(k)  *  arg  U(k)  +^(ln2  -  1),  (79) 

where 

arg  V(k)  =  -  J  ^  (arc  sin  ^  ,  (80) 


arg  U(k)  =  -  ^2  (^rc  sin-^-  .  (81) 

Therefore 

arg  (k)  =  x(ln  2  -  1)  -  (-)*'  (arc  sin-^  -  ^  .  (82) 

The  last  sum  can  be  written  in  the  form  arc  sin(x/2)  -  arc  sin(x/3)  -  x/6  +R(x),  where 

RW.J^(-)"(arcsln2-^).  (83) 


The  reason  for  this  separation  is  that  R(x)  is  small.  In  fact,  while  x  varies  from  1.01  to 
1.35,  R(x)  varies  from  .002  to  .004  and  the  entire  expression  (82)  varies  from  .33  to  .47. 


R(x)  is  calculated  by  expanding  arc  sin  (x/n)  in  an  infinite  series  and  then  interchanging 
the  order  of  summations.  This  is  permissible  since  all  series  involved  are  absolutely 
convergent.  From 


.  X 


arc  sin— 
n 


(84) 


it  follows  that 


The  summations  can  now  be  carried  out  since  the  function  ti  (s)  =  S  (-)"  n“®  is  related 

nz4 

lo  the  Zeta  function  of  Riemann  and  its  values  for  s  =  3,  5,  7,  etc.  can  be  found  with  the  aid 
of  the  tables  of  the  latter. 
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Table  of  the  Auxiliary  Function 
i1(s) 

s  3  5  7  9  11 

T)  .010  4943  .000  7454  .000  0504  .000  0034  .000  0002 

The  numerical  form  of  R(x)  is: 

R(x)  =  (17490x*  +  559x®+  22x’  +  x®+  ...)  •  10*’  ,  (86) 

while 

arg  K^(k)  =  -  .140  1861  x  -  arc  sin -^  +  arc  sin  R(x).  (87) 

In  order  to  achieve  the  necessary  symmetry  for  the  calculation  of  p"  for  norma] 
incidence,  the  origin  has  to  be  moved  to  a  point  half-way  between  the  edges,  such  as  the 
point  Q  in  Figure  2.  When  this  is  done  the  expression  for  the  phase  of  transmission 
coefficient  becomes  in  place  of  (21): 

t'  =  arg  K^{k  cos  ff)  -  arg  K^(k).  (88) 

Therefore,  in  view  of  (59)  and  (76) 

p«  =  2-7^  -  p'  +  ff  »  -2K^(x)  +  ff.  (89) 

It  is  therefore  necessary  to  calculate  K+(k).  For  normal  incidence  the  procedure  is 
similar  to  the  calculation  of  Klf(if).  However,  in  the  present  case  the  convenient  auxiUary 
variable  is  y  =  nx  =  KaA  . 

From: 


Hence 


arg  K^(k)  =  y(ln2  -  1)  (arc  sin  "-j)  • 


(92) 
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The  sum  is  again  split  as  follows: 


where 


arc  sin  *-  arc  sin  y/VF-  y/6  +  S(y), 


S(y)  =  S  (-)”  (arc  sin 
'''  n=4  '  \  i/n*  -  1 


“  i  ('  -T- 

The  binomial  theorem  gives: 

*  y«/,  5  -2  35  -4  \ 

These  expansions  may  now  be  substituted  in  the  series  of  arc  sin  u.  Absolute  conver¬ 
gence  of  all  series  permits  a  rearrangement  of  the  terms  resulting  in  the  following  expression: 

+  y®  -fn*®  +-|  n'®  +...|+  (96) 


In  a  compact  form 


^ere 


40  „?4  2 

S(y)  =  g^y  +  |g3y»  +  ^  g,y*+.... 


gj  =  I TJ  (3)  +  I  n  (5)  +  ^  ?7  (7)  +...  =  .005  5434, 
g,  =  0  (3)  +-§  (5)  +  ^ »]  (7)  +...  =  .011  7153. 

Similar  expressions  give  for  ggy  g,  and  g^,  8878,  679  and  46  ten-millicmths  respectively. 
The  numerical  form  of  the  remainder  is: 

S(y)  =  (55434y  +  19525y*  +  666y®  +  30y’  +  2y*  +  ...)•  10"’. 


(98) 


26  NAVAL  RESEARCH  LABORATORY 

\ 

The  final  expression  for  K+(tc)  is 

arg  K^(k)  =  -arc  sin  y/Ys  +  arc  sin  y/  /S  -.140  1861  -S(y).  (99j 

In  order  to  obtain  p"  for  6^0,  values  of  K+{tc)  for  oblique  incidence  are  required. 
These  were  obtained  by  calculating  the  difference  between  K+{ic,e)  and  K+(>f ,  0)  where  th* 
latter  is  the  function  designated  heretofore  by  K^.(x).  ^ 

Since  U(w)  does  not  depend  on  e 


arg  K^(k,  e)  -  arg  K^(ic,  0)  =  arg  V(ic,  6)  +  arg  W(»c,  0)  -  2  arg  V(k,  0).  (iqqj 
From  (13)  it  follows  that 


/.  ak  sin  oY  1 
V  ~  2im  I  ~  4n*  * 


(101) 


Let  V  =  (ak/ir )  sin  0  =  x  sin  0  ,  then 

arg  V(*,0)  =  -  E  (arc  sin  .  ^  . 

V(2n  -  v)"  -  1 

and 

arg  W(«.  6)  -  -  f  (arc  aln  _  , 

Introduce  the  function  fii(v)  =  arc  sin  y  |(2n  +  v)*  -  ij  The 

-E  f_(v)  +  f  (-v)  -  2f  (0)- 
nsi  n  n  n 


i)  • 

-i)  • 

right  side  of  (100)  then  becomes 


The  summand  is  the  second  difference  of  the  function  fn(v)*  K  is  equal  to  fn  (0)v* 
plus  terms  of  higher  order  in  v*.  But  fn  (0)  is  approximately  y/4n* ,  when  1  and  y*  can 
be  neglected  in  comparison  to  4n® .  In  the  range  of  interest  y<l.  Therefore,  when  N  is  such 
that  1  «4N® ,  the  following  approximation  holds 


n:N  n 


(104) 


The  terms  corresponding  to  n  =  1,  2,  3  and  4  were  computed  on  the  automatic  comput¬ 
ing  machine  of  the  Operational  Research  Branch  for  angles  of  incidence  0  ranging  from  0 
to  35®  in  5®  intervals.  The  remainder  term 

^  Ej  n-®  =  .006099yv®  (105) 


was  then  added.  The  values  of  p"  -  tt  =  -2  K+(x,0  )  are  tabulated  to  four  places,  except 
for  0  =  35®  in  which  case  only  three  place  accuracy  is  available. 


THE  TABLES  AND  THEIR  USE 

Table  1  contains  p  =  Ir  I  as  a  function  of  x  =  2a/x  from  1.01  to  1,35  and  the  angle  of 
incidence  0  from  0®  to  60®.  In  Table  2  the  values  of  p"-  jr  are  tabulated  for  x  =  1.01 
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to  1.35  and' for  ^  =  tf  to  35°,  the  approximations  adopted  not  being  satisfactory  for  large 
values  of  9 . 

A  more  extensive  table  was  prepared  for  the  case  of  normal  incidence  with  values  of 
X  ranging  from  1.01  to  1.35.  This  presentation  (Table  3)  also  contains  the  functions  p'/2ir 
and  p"/2v  which  are  useful  in  practical  computations. 

As  an  example  for  the  use  of  these  functions  one  might  consider  the  computation  of 
the  phase  change  in  transmission  through  a  sheet. 

Equation  (72)  can  be  put  in  the  following  form; 


T'  o'  1  d  cos  6  \l/  ,  4/'  -  i}/ 


(106) 


Here 


4/  p"  nd 
2v  2v  '  ’ 


(107) 


since  the  rays  are  constrained  in  the  metal-plate  medium  so  that  Oi  =  90°.  Having  the 
angle  in  the  form  4//2ir  is  convenient,  for  4/  is  often  quite  large. 


The  difference  4/'  -  4/ is  generally  small.  It  can  be  computed  quickly  as  follows:  Let 
F  =  (1  +p?  )/(l  -  p*),  then  according  to  (71)  tan  4/'-T  tani^,  therefore 


4/'  -  t/>»tan(i/^'  -</')  = 


tan  4^'  -  tan  _ F  -  1 _ 

1  +  tan  p'  t3tn4'~  cotan  4/  +  F  tan  4^ 


(108) 


provided  F  -  1  is  small. 

No  detailed  tables  were  prepared  for  the  power  division  at  the  air -metal -plate  inter 
face.  Figure  5  contains  the  available  information. 


Table  4  gives  the  direction  of  the  first  diffracted  beam  as  a  function  of  x  and  the  angle 
of  incidence.  It  also  contains  the  angle  Oq  =  arc  cos  ti  as  a  fimction  of  x.  This  is  the  an¬ 
gle  of  incidence  for  which  the  interesting  equidistribution  of  power  noticed  by  L.  J.  Chu 
takes  place. 
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TABLE  1 

Magnitude  of  the  Reflection  Coefficient  p  for  Varying  Angle  of  Incidence 


x  =  -2a 

o 

It 

5° 

O 

O 

— 

15“ 

O 

O 

25“ 

30P 

1,01 

.7538 

.7530 

.7505 

.7462 

.7401 

.7318 

.751  ft 

1.02 

.6708 

.6697 

.6665 

.6611 

.6533 

.6428 

.R5qo 

1.03 

.6134 

.6122 

.6086 

.6025 

.5937 

.5818 

■5666 

1.04 

.5690 

.5677 

.5638 

.5572 

.5476 

.5348 

1.05 

.5327 

.5313 

.5272 

.5201 

.5100 

.4965 

.4792 

1.06 

.5019 

.5005 

.4961 

.4888 

.4782 

.4642 

.4462 

1.07 

.4752 

.4737 

.4692 

.4617 

.4508 

.4362 

.4176 

1.08 

.4517 

.4502 

.4456 

.4378 

.4266 

.4117 

.3926 

1.09 

.4307 

.4292 

.4245 

.4165 

.4050 

.3898 

.3704 

1,10 

.4118 

.4102 

.4055 

.3973 

.3857 

.3702 

.3504 

1.11 

.3947 

.3931 

.3882 

.3799 

.3681 

.3524 

.3323 

1.12 

.3790 

.3774 

.3724 

.3640 

.3520 

.3361 

.3158 

1.13 

.3646 

.3629 

.3579 

.3494 

.3373 

.3212 

.3006 

1.14 

.3512 

.3495 

.3445 

.3359 

,3237 

,3074 

.2867 

1.15 

.3388 

.3371 

.3321 

.3234 

.3110 

.2946 

.2737 

1.16 

.3273 

.3256 

.3205 

-3118 

.2993 

,2827 

.2617 

1.17 

.3166 

.3148 

.3096 

.3009 

.2883 

.2716 

.2504 

1.18 

.3065 

.3047 

.2995 

.2907 

.2780 

.2612 

.2399 

1.19 

.2970 

.2952 

.2900 

.2811 

.2684 

.2515 

.2301 

1.20 

.2880 

.2863 

.2810 

.2721 

.2592 

.2423 

.2208 

1.21 

.2796 

.2778 

.2725 

.2635 

.2507 

.2336 

.2120 

1.22 

.2716 

.2698 

.2645 

.2555 

.2426 

.2255 

.2038 

1.23 

.2640 

.2622 

.2569 

.2478 

.2349 

.2177 

.1959 

1.24 

.2568 

.2550 

.2497 

.2406 

.2276 

.2103 

.1885 

1.25 

.2500 

.2482 

.2428 

.2337 

.2206 

.2034 

.1815 

1.26 

.2435 

.2417 

.2363 

.2271 

.2140 

.1967 

.1748 

1.27 

,2373 

.2355 

.2300 

.2209 

.2077 

.1904 

.1684 

1.28 

.2314 

.2295 

.2241 

.2149 

.2017 

.1843 

.1623 

1,29 

.2257 

.2239 

.2184 

.2092 

.1960 

.1785 

.1564 

1.30 

.2203 

.2185 

.2130 

.2037 

.1905 

.1730 

.1509 

1.31 

.2151 

.2133 

.2078 

.1985 

.1852 

.1677 

.1455 

1.32 

.2101 

.2083 

.2028 

.1935 

.1802 

.1626 

.1404 

1.33 

.2053 

.2035 

.1980 

.1887 

.1754 

.1578 

.1355 

1.34 

.2007 

.1989 

.1934 

.1841 

.1707 

.1531 

♦.1737 

1.35 

.1963 

.1945 

.1890 

.1796 

.1662 

.1486 

♦.1980 

♦Second  lobe  present. 
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TABLE  1 
-continued- 


o 

O 

45° 

.6903 

.6687 

.5908 

.5641 

.5235 

.4938 

.4722 

.4405 

.4306 

.3974 

.3957 

.3614 

.3658 

.3306 

.3395 

.3037 

.3163 

.2798 

.2955 

.2585 

.2767 

.2393 

.2596 

.2218 

.  .2439 

.2059 

.2294 

.1912 

.2161 

.1776 

.2037 

.1650 

.1921 

.1533 

.1813 

♦.2227 

.1712 

♦.2582 

.1617 

♦.2845 

.1528 

♦.3062 

♦.1793 

♦.3248 

♦.2215 

♦.3412 

♦.2481 

♦.3559 

♦.2691 

♦.3693 

♦.2870 

♦.3815 

♦.3026 

♦.3928 

♦.3165 

♦.4033 

♦.3292 

♦.4131 

♦.3408 

♦.4223 

♦.3515 

♦.4310 

♦.3615 

♦.4392 

♦.3709 

♦.4469 

♦.3797 

♦.4543 

♦.3880 

♦.4614 

50° 

55° 

.6415 

.6068 

.5307 

.4886 

.4569 

.4107 

.4012 

.3524 

.3565 

.3058 

.3193 

.2672 

.2875 

.2344 

.2598 

.2059 

.2353 

.1808 

.2135 

♦.1849 

.1939 

.1761 

.1598 

♦.3479 

♦.2347 

♦.3747 

*.2776 

♦.3973 

♦.3085 

*.4169 

♦.3335 

♦.4342 

♦.3547 

♦.4497 

♦.3732 

♦.4638 

*.3896 

♦.4766 

.7074 

.6122 

.5474 

.4978 

.4575 

.4236 

.3944 

.3688 

.3461 

.3258 

.3073 

.2905 

.2751 

.2609 

.2478 

.2356 

.2242 

.2135 

.2036 

.1942 

.1853 

.1770 

.1690 

.1615 

.1544 

.1476 

.1412 

♦.1948 

*.2206 

♦.2402 

♦.2566 

♦.2709 

♦.2837 

♦.2953 

♦.3060 


♦.4621 

♦.4713 

♦.4800 


*.5033 

♦.5104 

♦.5171 

♦.5234 

*.5296 


♦.5187 

♦.5275 


♦.5581 

♦.5647 

♦.5711 

♦.5772 


*®®cond  lobe  present 


60° 

.5616 

.4346 

.3521 

.2909 

.2423 

.2024 

.1686 

♦.2801 

♦.3370 

♦.3771 

♦.4086 

♦.4345 

♦.4566 

♦.4758 

♦.4927 

♦.5078 

♦.5214 

♦.5338 

♦.5452 

♦.5557 

♦.5654 

♦.5744 

♦.5829 

♦.5908 

♦.5982 

♦.6053 

♦.6120 

♦.6183 

♦.6244 

♦.6300 

♦.6356 

♦.6408 

♦.6458 

♦.6506 

♦.6553 
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TABLE  2 

Phase  of  the  Back  Surface  Reflection  Coefficient 


Table  of  P"  -  <r  (in  radians) 


X 


1.04 

1.05 

1.06 

1.07 

1.08 

1.09 

1.10 

1.11 

1.12 

1.13 

1.14 

1.15 

1.16 

1.17 

1.18 

1.19 

1.20 

1.21 

1.22 

1.23 

1.24 

1.25 

1.26 

1.27 

1.28 

1.29 

1.30 


6  =  0“ 

n 

o 

o 

15“ 

20“ 

25“ 

.1049 

.1064 

.1107 

.1180 

.1284 

.1425 

.1490 

.1511 

.1573 

.1679 

.1832 

,2038 

.1832 

.1858 

.1937 

.2070 

.2265 

.2527 

.2124 

.2155 

.2248 

.2408 

.2639 

.2952 

.2384 

.2419 

.2527 

.2710 

.2977 

.3340 

.2622 

.2662 

.2782 

.2990 

.3292 

.3704 

.2843 

.2888 

.3023 

.3253 

.3590 

.4051 

.3052 

.3101 

.3249 

.3503 

.3875 

.4586 

.3250 

.3303 

.3465 

.3742 

.4150 

.4712 

.3440 

.3498 

.3673 

.3975 

.4418 

.5032 

.3623 

.3685 

.3874 

.4200 

.4681 

.5348 

.3799 

.3866 

.4069 

.4420 

.4939 

.5662 

.3970 

.4042 

.4259 

.4636 

.5194 

.5976 

.4137 

.4213 

.4445 

.4848 

.5447 

.6290 

.4300 

.4381 

.4628 

.5058 

.5700 

.6605 

.4460 

.4546 

.4809 

.5265 

.5951 

.6923 

.4616 

.4708 

.4986 

.5472 

.6202 

.7245 

.4770 

.4866 

.5162 

.5677 

.6454 

.7571 

.4921 

.5022 

.5335 

.5881 

.6707 

.7902 

.5071 

.5177 

.5508 

.6085 

.6964 

.8240 

.5218 

.5331 

.5678 

.6289 

.7222 

.8587 

.5364 

.5482 

.5849 

.6493 

.7483 

.8941 

.5508 

.5632 

.6018 

.6698 

.7747 

.9305 

.5650 

.5782 

.6187 

.6904 

.8015 

.9680 

.5792 

.5930 

.6356 

.7113 

.8288 

1.0069 

.5932 

.6077 

.6525 

.7323 

.8567 

1.0471 

.6072 

.6223 

.6693 

.7532 

.8850 

1.0888 

.6211 

.6369 

.6862 

.7744 

.9140 

1.1323 

.6349 

.6515 

.7031 

.7960 

.9437 

1.1781 

.6486 

.6660 

.7200 

.8178 

.9743 

1.2261 
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TABLE  3 

Tables  for  Normal  Incidence 


*  X 

n 

P 

P' 

lit 

P“-  IT 

P" 

lit 

1.01 

.14037 

.7538 

-.6589 

-.10487 

.1049 

.51670 

1.02 

.19706 

.6708 

.1490 

.52371 

1.03 

.23959 

.6134 

-.6738 

-.10724 

.1832 

.52916 

1.04 

.27467 

.5690 

-.6813 

-.10844 

.2124 

.53380 

1.05 

.30491 

.5327 

-.6889 

-.10964 

.2384 

.53794 

1.06 

.33167 

.5019 

-.6955 

-.11068 

.2622 

.54173 

1.07 

.35575 

.4752 

-.7041 

-.11206 

.2843 

.54525 

1.08 

.37770 

.4517 

-.7117 

-.11327 

.3052 

.54857 

1.09 

.39789 

.4307 

-.7194 

.3250 

.55173 

1.10 

.41660 

.4118 

-.7271 

-.11573 

.3440 

.55475 

1.11 

.43403 

.3947 

-.7349 

.3623 

.55766 

1.12 

.45034 

.3790 

-.7427 

.3799 

.56047 

1.13 

.46567 

.3646 

-.7506 

.3970 

.56319 

1.14 

.48014 

.3512 

-.7585 

.4137 

.56585 

1.15 

.49382 

.3388 

-.7664 

.4300 

.56844 

1.16 

.50679 

.3273 

-.7744 

-.12324 

.4460 

.57098 

1.17 

.51912 

.3166 

-.7824 

-.12452 

.4616 

.57347 

1.18 

.53086 

.306B 

-.7905 

-.12581 

.57592 

1.19 

.54207 

.2970 

-;7986 

-.12710 

.57832 

1.20 

.55277 

.2880 

-.8067 

-.12840 

.5071 

.58070 

1.21 

.56302 

.2796 

-.8150 

-.12971 

.5218 

.58305 

1.22 

.57283 

.2716 

-.8232 

-.13102 

.5364 

.58536 

1.23 

.58225 

.2640 

-.8316 

-.13235 

.5508 

.58766 

1.24 

.59130 

.2568 

-.8399 

-.13368 

.5650 

.58993 

1.25 

.60000 

.2500 

.  -.8484 

-.13502 

.5792 

.59218 

1.26 

.60837 

.2435 

-.8569 

-.13637 

,  .5932 

.59442 

1.27 

.61644 

.2373 

-.8654 

-.13773 

.6072 

.59664 

1.28 

.62422 

.2314 

-.8740 

-.13910 

.6211 

.59884 

1.29 

.63172 

.2257 

-.8827 

-.14048 

.6349 

.60104 

1.30 

.63897 

.2203 

-.8914 

-.14187 

.6486 

.60323 

1.31 

.64597 

.2151 

-.9002 

-.14327 

.6623 

.60540 

1.32 

.65275 

.2101 

-.9091 

-.14468 

.6759 

.60758 

1.33 

.65930 

.2053 

-.9180 

-.14610 

.6895 

.60974 

1.34 

.66564 

.2007 

-.9270 

-.14753 

.7031 

.61191 

1.35 

.67179 

.1963 

-.9360 

-.14898 

,7167 

.61406 
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TABLE  4 

Azimuth  of  the  First  Diffracted  Beam  and  the  Angle 
=  arc  cos  b. 


Values  of  9'  in  degrees 


x=-^ 

*  X 

^0 

V  Alt 

5^1 CCO  ^ 

9  =  35“ 

9  =  40“ 

^  =  45“ 

9  =  50“ 

e  =  55“ 

9  =  60“ 

9\ 65* 

1,05 

1.06 

72.25 

70.63 

86.81 

1.07 

69.16 

78.66 

1.08 

67.81 

80.34 

74.33 

71  f\A 

1.09 

66.55 

75.66 

68.21 

1.10 

65.38 

87.48 

72.20 

65.77 

1.11 

64.28 

79.31 

69.35 

63.57 

1.12 

63.23 

75.14 

66.88 

61.57 

1.13 

62.25 

71.95 

64.67 

59.72 

1.14 

61.31 

81.24 

69.27 

62.67 

58.00 

1.15 

60.41 

76.68 

66.92 

60.82 

56.39 

1.16 

59.55 

73.35 

64.82 

59.11 

54.87 

1.17 

58.73 

70.62 

62.91 

57.50 

53.43 

1.18 

57.94 

81.04 

68.26 

61.14 

55.99 

52.03 

1.19 

57.18 

76.80 

66.15 

59.49 

54.55 

50.75 

1.20 

56.44 

73.65 

64.24 

57.94 

53.19 

49.50 

1.21 

55.74 

71.05 

62.48 

56.49 

51.89 

48.30 

1.22 

55.05 

85.24 

68.79 

60.84 

55.10 

50.65 

47.14 

1.23 

54.39 

79.49 

66.77 

59.31 

53.79 

49.46 

46.03 

1,24 

53.75 

75.96 

64.93 

57.87 

52.54 

48.32 

44.96 

1.25 

53.13 

73.18 

63.24 

56.51 

51.34 

47.22 

43.92 

1.26 

52.53 

70.82 

61.67 

55.21 

50.19 

46.16 

42.92 

1.27 

51.94 

68.75 

60.19 

53.98 

49.08 

45.14 

41.95 

1.28 

51.38 

81.47 

66.88 

58.80 

52.79 

48.03 

44.15 

41.01 

1.29 

50.82 

77.64* 

65.18 

57.49 

51.66 

46.99 

43.19 

40.10 

1.30 

50.28 

74.77 

63.60 

56.24 

50.57 

46.00 

42.26 

39.21 

1.31 

49.76 

72.39 

62.12 

55.05 

49.52 

45.04 

41.35 

38.35 

1.32 

49.25 

70.32 

60.73 

53.91 

48.51 

44.11 

40.48 

37.51 

1.33 

48.75 

68.46 

59.43 

52.81 

47.54 

43.20 

39.62 

36.62 

1.34 

48.27 

66.77 

58.18 

51.76 

46.59 

42.33 

38.79 

35.89 

1.35 

47.79 

65,22 

57.00 

50.75 

45.68 

41.48 

37.99 

35.11 
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APPENDIX 

Calculation  of  the  Magnitude  of  the  Reflection  and  Transmission  Coefficients 


In  order  to  find  the  magnitudes  of  certain  coefficients  occurring  on  pages  8-11,  it  is 
necessary  to  evaluate  the  absolute  value  of  the  complex  function  K^.(w)  and  its  derivative 
for  certain  special  values  of  the  variable  w.  In  particular,  |K^(w)  i  is  needed  for  w=kcos0 
and  K,  while  lK+(w)  I  is  needed  for  w  =  k  cos  6  and  -k  cos  9' .  Direct  substitution  in  the 
infinite  products  defining  K^(w)  leads  to  extremely  cumbersome  expressions*,  a  somewhat 
different  method  is  therefore  preferable. 


The  functions  K^{w)  and  K_(w)  are  defined  by  (9)  and  (15).  On  multiplying  their 
absolute  values  the  following  equation  results: 


a 

""  2 

w  +  k  cos  9  w  -  tc 

V(w) 

m 

,w)U(-w)| 

w  -  k  cos  0  w  +  K 

V(-w)| 

lW(-w)U{w)  1 

K^(w)K.(w)J  =1 

In  view  of  (8)  the  left  side  of  (lA)  might  also  be  written  as  |f(w)j^K^(w)]  *  L 


(lA) 


When  w  is  real,  the  last  factor  on  the  right  of  (lA)  is  always  one,  since  W(-w)  =  W(w)* 
and  U(-w)  =  U(w)*.  When  is  real  then  also  V(-w)  =  V(w)*  and  in  this  case 


If(w)|  lK^(wf  =  I 


w  +  k  cos  9  w  -K 
w  -  k  cos  0  w  +»t  * 


When  Aj  is  imaginary  but  A j  is  real,  for  every  real  w 


(2A) 


and  therefore 


V(wl 

k  cos  0 '  +  w 

V(-w) 

k  cos  0 '  -  w 

Irr  ^ 

w  +  k  cos  0]  jw  -  K 

w  +  k  cos  0' 

2 

w  -  k  cos  01  Iw  +  K 

w  -  k  cos  0’ 

(3A) 


To  derive  an  expression  for  tK^(k  cos  0)|,  it  is  only  necessary  to  multiply  (2A)  or  (3A) 
*>y  (w  -  k  cos  9)  and  then  let  w  tend  to  k  cos  9.  According  to  L'  Hospital’s  rule 


lim  f(w)(w  -  k  cos  9)  =  l/ak  cos  9, 

w-4  k  cot 9 

^^refore 

Ik  (k  cos  0)1*=  a*k*cos*  0 
+  '  k  cos  9  +  K 

^hen  A^  is  real  and 


lK^(k  cos  0)1*=  a*k*  cos*  0 


k  cos  0-K  cos  0  +  cos  9' 
k  cos  0+K  cos  0  -  cos  0' 


(4A) 


(5A) 


^hen  Aj  is  imaginary. 
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The  omission  of  the  absolute  value  signs  on  the  right  of  (4A)  and  (5A)  is  justified  in  vi 
of  the  fact  that  cos  0  >  n  for  real  Aj  and  that  for  imaginary  Ai,  n  =  cos  lies  between*^ 
cos  6  and  cos  0'  as  seen  from  (41).  “ 


The  calculation  of  [K^(*f)j  requires  the  evaluation  of  the  limit; 


sin  a  Vk*  -  w“ 


-  A®  cos  a  Vk®  -A®  -  cos  (ak  sin  9)  *■** 


Since  k®  -  a  ®  =  ff®/a®  this  limit  becomes 


1  +  cos(ak  sin  0) 


Therefore 


|k^(*)  I*  =  k  cos  0  -K  [l  +  cos  slo  0)] 


when  Aj  ia'  real  .and 


It,  /  a  h  JT®  k  cos  0  +A  k  cos  0'+K  c.  . -i 

|K^(*)  I  =  -  4^  k  cos  e'-tc  e  +  ‘=0® 

when  Aj  is  imaginary. 

Since  K^.(w)  has  a  simple  root  at  w  =  -  k  cos  0 

.  K  (w) 

|K'(-k  cos  0)  s  lim  - ± ^  . 

'  +  '  w^kco*e  w  +  kcos0 


Equations  (2A)  and  (3A)  are  divided  by  (w  +  k  cos  0)  and  the  left  side  is  written  in  the 
form 

r  K  (w)  1* 

«*)(W.kC05«[— 


From  {w  +  k  cos  fl)  f(w)  =  - 1 A  k  cos  e  then  follows: 

cos^ 

|K'(-k  cos  0)1®  =  I-* 

'  +  ''4  k  cos  0  -A 

when  A  ^  is  real  and 


|K'(-k  cos  0)1®  =  a®^.cose.^  cos  A  ^  cos  .0' 
+  ''4  k  cos  0-A  cos  0  +  cos  9 

when  A  ^  is  imaginary. 

The  calculation  of  K^(-k  cos  0')  is  similar.  In  this  case 
lim  f(w)  (w  +  k  cos  0')  =  -  1/a  k  cos  0'. 

w^Thcos  9' 
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Hence 


Ir,  /  1  «/\l2  a*  k  cos  6'+  K  cos  6'  -  cos  9 

lK^{-k  cos  0  )|  -  4  gQg  g/_  K  cos  6'  +  cos  6 

From  (23),  (4A)  and  (9A)  follows 


in 


k  cos  6  -K 


k  cos  B  +K 

for  real  Aj,  while  for  imaginary  (23),  (5A)  and  (lOA)  give 

k  cos  0  -  A  cos  0  +  cos  0' 


in  = 


k  cos  0  +A  cos  0  -  cos  0' 


(llA) 


(12A) 


(13A) 


Expressions  for  it!  and  r^  are  obtained  in  the  same  way.  The  formulas  (20),  (37)  and 
(39)  result. 


♦  ♦  * 


! 
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